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Abstract. This paper is devoted the the study of the mean field limit 
for many-particle systems undergoing jump, drift or diffusion processes, 
as well as combinations of them. The main results are quantitative es- 
timates on the decay of fluctuations around the deterministic limit and 
of correlations between particles, as the number of particles goes to in- 
finity. To this end we introduce a general functional framework which 
reduces this question to the one of proving a purely functional estimate 
on some abstract generator operators [consistency estimate) together 
with fine stability estimates on the flow of the limiting nonlinear equa- 
tion [stability estimates). Then we apply this method to a Boltzmann 
collision jump process (for Maxwell molecules), to a McKean-Vlasov 
drift-diffusion process and to an inelastic Boltzmann collision jump pro- 
cess with (stochastic) thermal bath. To our knowledge, our approach 
yields the first such quantitative results for a combination of jump and 
diffusion processes. 
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1. Introduction 



Fundamental in Boltzmann's deduction of the equation bearing his name is 
the " stosszahlansatz", or chaos assumption. Vaguely expressed this assump- 
tion means that when two particles collide, they are statistically uncorrelated 
just before the collision. After the collision they are not, of course, because, 
for example, the knowledge of the position and velocity of one particle that 
just collided gives some information on the position of the collision partner. 
And while the correlations created by collisions decrease with time, they 
never vanish in a system of finitely many particles, and hence the Boltz- 
mann assumption could only be true in the limit of infinitely many particles. 

A mathematical framework for studying this limit is the so-called BBGKY 
hierarchy (see Grad [17], Cercignani [9] and the book |11) by Cercignani et 
al), which consists of a family of Liouville equations, each describing the 
evolution of an A^-particle system (deterministic, in this case), and whose 
solutions are densities in the space of A^-particle configurations in phase 
space. The BBGKY hierarchy describes in a systematic way the evolution 
of marginal distributions. Formally, and under appropriate assumptions, 
most notably the chaos assumption, the one-particle marginal of solutions 
to the A^-particle Liouville equation, converges to solutions of the Boltzmann 
equation. 

It would take until 1974 before a mathematically rigorous proof of this 
statement was given by Lanford |27) . While this is a remarkable result, it 
only proves that the Boltzmann equation is a limit of the A'^-particle systems 
for a fraction of the mean free time between collisions, and this is essentially 
where the problem stands today (see however [22] for a large time limit in 
a near to the vacuum framework; it is worth emphasizing that in such a 
framework no more collisions occur than in Lanford's framework). 

In order to avoid some of the difficulties related to the deterministic evo- 
lution of a real particle system, Kac [23] invented a Markov process for a 
particular A^-particle system, and gave a mathematically rigorous definition 
of propagation of chaos. He then proved that this holds for his Markov sys- 
tem, and thus obtained a mathematically rigorous derivation of a simplified 
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(spatially homogeneous) Boltzmann equation in this case, usually called the 
Kac equation. 

Kac's work provides the framework of this paper, and we will now describe 
our main results. We let E be the state space of one particle (usually W^, 
but metric, separable and locally compact is fine). A sequence of probability 
measures (/^)jv=o> where each G V{E^) is symmetric in the sense that 
it is invariant under permutation of the coordinates, is said to be f -chaotic 
for some probability measure / G V{E) if for each A; > 1 and functions 
(pj G Ci,{E),j = 1, . . . ,k, (continuous bounded), 

n k k n 

(1.1) hm / r[0,(z,)/^(dzi,...,dzjv) = T\ <l),{z)f{dz). 

We next consider a family of time dependent probability measures (//^)^=0' 
being the distributions of the states of Markov processes in E^ . The Markov 
process is said to propagate chaos if given an initial family of A/^-particle 
distributions (/j>[);y=0' ^^^^ /m-chaotic, there is a time dependent distri- 
bution ft such that (//^)^=x is /^-chaotic. In this paper we are interested 
in specific equations that govern the evolution of ft] but it is important to 
bear in mind that they are in general nonlinear, and it may be difficult to 
prove well-posedness in function spaces relevant for proving the propagation 
of chaos. 

The main results of this paper are abstract. We consider: 

• the family of A^-particle systems represented by a Markov processes 
{Z^)t>o in some product space E^ , with = {Zi^It ■ ■ , ^N,t), 
and the corresponding probability distributions {ft^)^^i, solving 
Kolmogorov's backward equation^, and where the distributions //^ 
belong to a suitable subspace of V{E^), 

• a (nonlinear) equation defined on a subspace of V{E), which is the 
formal limit of the equations governing one-particle marginals of f^: 

(1.2) ^/t = Q(/i), k.eV{E). 
Then we: 

• provide conditions on the processes and related function spaces that 
guarantee that (//^)7v=i is /t-chaotic for t > 0, 

• give explicit estimates of the rate of convergence in (II. ID : more pre- 
cisely, for any T > 0, £ G N* and 4>j £ F C Cb{E) (j = 1, . . . ,£), 



As we will see in the following section, the formalism also works if Zt satisfies a 
deterministic evolution, in which case Liouville's equations replace Kolmogorov's backward 
equation 
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there is a constant e(A^) converging to zero as — )> oo such that 

(1.3) 



sup 

te[o,T] 



I f{Uz,)f!'{dzu...AzN)-\{l Hz)ft{dz) 

Je^ j=i j=i -'^ 



which holds for N >2£ and a suitably chosen space J^; if J-" is dense 
in Cb{E), this implies in particular the propagation of chaos. 

To this end, our starting point is a technique that goes back at least to 
Griinbaum |21| . which consists in representing an A^-particle configuration 
as a sum of Dirac measures, 



1 ^ 



and proving that, in a weak sense, |U^jv converges to //^. In fact, because 

is random, fi^j^ is a random measure in ViE), which has a probabil- 

ity distribution e V{V{E)). Proving the propagation of chaos is here 
equivalent to proving that ~^ ^ft 'Pi'^i^)) when — )• oo. The error 
e(A^) is dominated by, on the one hand, how well the initial measure /;„ can 
be approximated by a sum of N Dirac measures, and, on the other hand, 
estimates comparing the equations for //^ and ft- These estimates depend 
on rather technical assumptions, and although the abstract main theorem is 
stated in Section [21 the assumptions are stated in full detail only in Section 

With the main theorem of this paper in hand, proving the propagation of 
chaos for a particular A'^-particle system is reduced to proving: 

(i) a purely functional estimate on the dual generator of the A^- 
particle dynamics which establishes and quantifies that, at first or- 
der, is linked to the mean field limit generator Q {consistency 
estimate); 

(ii) some fine stability estimates on the flow of the mean field limit equa- 
tion involving the differential of the semigroup with respect to the 
initial data (stability estimates) . 

Point (i) of our method is largely inspired from the "duality viewpoint" 
of Griinbaum's paper |21j where he considered the propagation of the chaos 
issue for the Boltzmann equation associated to hard-spheres (unbounded) 
kernel. As he confessed himself the proof in [2Tj was incomplete due to the 
lake of suitable stability estimates, i.e. precisely the point (ii) of our method. 

It is worth emphasizing that after we had finished writing our paper, we 
were told about the recent book [26j by Kolokoltsov and his series of pa- 
pers on nonlinear Markov processes and kinetic equations. These interesting 
works focus on fiuctuation estimates of LLN and CLT types in the general 
framework of nonlinear Markov processes, and in some sense they generalise 
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to several other kinetic models the Griinbaum's duality viewpoint (although 
Kolokoltsov seems to not be aware of that earlier work). However we were 
not able to extract from these works a full proof in the cases when the gen- 
erator is an unbounded operator and weak distances have to be used. While 
the comparison of generators for the many-particle and the limit semigroup 
present in both |21| and |26| is reminiscent of our work, we believe that the 
main novelty of the present paper is to achieve, for the first time, both the 
fine stabilities estimates in point (ii) and the consistency estimate in point 
(i) in appropriate spaces (with weak topologies), in such a way that they 
may be combined and they lead to the already mentioned propagation of 
chaos result with quantitative estimates. 

We illustrate the method by proving the propagation of chaos for three 
different well known examples: 

(a) We first consider the Boltzmann equation for Maxellian molecules 
with angular cutoff. For such a bounded kernel case the result is 
well-known since the pioneering works of Kac [231 [24] and McKean 
[29j (who prove the propagation of chaos without any rate) and from 
the works by Graham and Meleard |18| [T9l [201 130| (where the au- 
thors establish the propagation of chaos with optimal rate 0{1/N)). 
In these papers, the cornerstone of the proof is a combinatorial argu- 
ment applied to the equation on the law (Wild sum expansion) or to 
the stochastic flow (stochastic tree). These approaches are restricted 
to a constant (or at least bounded) collision rate. 

(b) The second example is the McKean- Vlasov model. For such a model 
again, propagation of chaos is well-known and has been extensively 
studied. One of the most popular and efficient approach to deal 
with this model is the so-called "coupling method" introduced in the 
1970's, which yields the optimal convergence rate C(l/\/iV) (note 
that the difference between these two optimal rates in (a) and (b) 
comes from the fact that they are not measured with the same dis- 
tance). We refer to the lecture notes |37i i30j as well as to the ref- 
erences therein for a detailed discussion of that method. We also 
refer to [5j and the references therein for recent developments on the 
subject. 

(c) The third example is a mixed collision- diffusion equation which arises 
from granular gas modeling. For such a model, it seems that both 
the "combinatorics method" and the "coupling method" fail while 
our present method is robust enough to apply and yield quantitative 
chaos estimates. Let us also emphasize that the BBGKY method 
and the nonlinear martingale method (see again [371 [30] or [H I31| ) 
may also apply but would give a propagation of chaos without any 
rate since they are based on compactness arguments. 

Let us emphasize that it is not difficult to write a uniform in time version 
of Theorem 2.1: in short, if the assumptions (Al) to (A5) are satisfied with 



6 



S. MISCHLER, C. MOUHOT, AND B. WENNBERG 



T = +00, then the conclusion of the main abstract Theorem 12.11 holds with 
T = +00 and the proof is unchanged. But such an abstract theorem does 
not readily apply to the examples (a), (b) and (c) discussed above. More 
precisely, it is indeed possible to prove quantitative uniform in time propaga- 
tion of chaos by our method (for the elastic Boltzmann model for instance), 
but the price to pay is a significant modification to the set of assumptions 
(Al) to (A5). This issue is addressed in our companion paper [32] where 
the abstract method is developed in a more general framework in order to (1) 
apply it to Boltzmann collision models associated to unbounded collisions 
rates, (2) develop a theory of uniform in time propagation of chaos estimates. 
We shall consider the question of uniform in time chaoticity estimates for 
the McKean-Vlasov equation in future works. 

The plan of the paper is as follows. In Section [21 we present the method 
in an abstract framework, by first setting up a functional framework that is 
appropriate for comparing the A^-particle dynamics with the limiting dynam- 
ics, and we establish the abstract quantitative propagation of chaos (Theo- 
rem [2?TI- The main steps of the proof are given as well, but these rely on 
some technical assumptions and lemmas which are postponed to Section 31 
The functional framework is developed with the necessary details in Sec- 
tion [3l where we also develop a differential calculus for functions on V{E), 
as needed for studying the nonlinear semigroup. In Section [5l we apply 
the method to the Boltzmann equation associated to the Maxwell molecules 
collision kernel with Grad's cut-off. In Section [6l we apply the method to 
the McKean-Vlasov equation, and finally, in Section [71 it is applied to some 
mixed jump and diffusion equations motivated by granular gases. 

2. Propagation of chaos for abstract TV-particle systems 

In this section we introduce the mathematical notation used in the paper, 
make precise statements of the results, and describe the main steps of the 
proof, leaving the details of the proofs to the next sections. 

2.1. The A^-particle system. The phase space of the A^-particle system 
isi E^/6n- Here E is assumed to be a locally compact, separable metric 
space, and &N denotes the symmetric group of order N. This means that 
we identify all points in the A^-particle phase space that can be obtained 
by permutation of the particles, so that if Z = (zi, . . . ,zn) G E^ we 
have (zi, . . . , zn) ~ i^ai, ■ ■ ■ , Za^), where (cti, . . . , ajy) is any permutation 
of {1, . . . , N}. The evolution in phase space may be a stochastic Markov 
process or the solution to an Hamiltonian system of equations. In both 
cases, we denote by (2^/^)j>q the flow of the process. 



The phase space of a reahstic A^'-particle system may be a subspace of the A'^-fold 
product, determined e.g. by energy constraints or by the fact that particles of finite size 
may not overlap. However, in the hmit — 00, all of i5 should be accessible for any 
given particle. 
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Figure 1. A summary of spaces and their relations. Semi- 
groups are in most cases given together with their generators, 



as in S^\A. 



Figure [T] illustrates the relation between the different objects that we con- 
sider. The A'^-particle system is represented in the upper left corner of Fig- 
ure [TJ The different mathematical objects in this diagram are explained 
along the following subsections. 

2.2. Master equations, Liouville's equations and their duals. Let 

'Psym{E^) denote the proability measures on that are invariant under 
permutation of the indices in Z = (zi, . . . , zn) G E^ . The flow induces 
a semigroup of operators 5^ on Vsym{E^) defined through the formula 

(2.1) V/^EP,y.„(i?^), (^Ga(ii;^), 

where the bracket denotes the duality bracket between V{E^) and Ci,{E^): 

= / HZ)f{dZ), 
Jen 

and ^Zin denotes the conditional expectation with respect to the initial con- 
dition Zl^ = Zi„. This semigroup is the solution to Kolmogorov's forward 
equation in the case where {Z^)t>o is a random process (this equation is 
often called the master equation) , and of the Liouville equation in the Hamil- 
tonian case. We always assume that S^^ preserves the symmetry under per- 
mutation, and therefore restricts to an evolution semigroup on Vsyin{E^). 
There is a dual semigroup of Sf^ , that acts on Cb{E^), the set of bounded 
continuous functions on . We write this semigroup T/^, and denote its 
generator by . The two semigroups are related by 

(2.2) (/^^,^0) = (5f/^<A>. 
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Markov processes such that (1) is a contraction in Cq{E ), the set of 
continuous functions that vanish at infinity, and (2) 1 1— )■ T^^cj) is continuous 
for any (/) G Cq{E^), are known as Feller processes. 

Hence the upper part of the diagram represents a (random) process, and 
its distributions. The lower part of the diagram essentially shows the same 
thing as induced by the map /i^, as we shall now see. 



2.3. The limiting dynamics. The components (zi, . . . , zn) of Z e E /G^ 

represent the positions (in generalized sense, i.e. in thephase space E) of 
the N particles. These N particles can also be uniqueljo represented as an 
empirical measure, that is a sum of Dirac measures: 

1 ^ 

(2.3) Z = (zi,...,zjv) ^ /^z=^E^-r 

The resulting measure is normalized so as to give a probability measure, 
which is obviously independent of any permutation of the indices. The set 
of such empirical measures is denoted by Vn{E)- Probability measures on 
E are denoted by V{E), so that VNiE) C ViE). 

When the number of particles go to infinity, we may have /i^ — ?■ / G V{E), 
where now / is a distribution of particles in E. We call a "limiting equation 
for the iV-particle systems" the (usually nonlinear) equation of the form (II. 2p 
that is satisfied by the probability distribution ft obtained as the limit of 
fi^ , and we write its solution in the form of a nonlinear semigroup S^^: 

NLi 



is the solution of 



ft = (/in) 



dtft = Qift) , /O = /in • 



The main result of this paper can be seen as a perturbation result: given a 
solution to the limiting equation, we consider a sequence of measures IJ,2. G 
Vn{E) such that 

and prove that for all t in some interval < t < T, 

The convergence is established in the weak topology for the law of the random 
empirical measures, as will be explained next. 



'The equivalence of particle configurations under permutation is used here. 



NEW APPROACH TO QUANTITATIVE PROPAGATION OF CHAOS 



9 



2.4. A^-particle dynamics of random measures and weak solutions of 
the limiting equation. A random point Z G with law G VaymiE^) 
can be identified with a random measure, denoted by G Vn{E), whose 
law is induced from . We denote this law TTpf^ G V{V{E)). Note that 
since ii^ is a separable metric space, then so is V{E) by Prokhorov's Theorem, 
and we may define the space V{V{E)) of probability measures on V{E), as 
well as the set of continuous bounded functions, denoted Ch{V{E)), which is 
the dual of V{V{E)). In Section [3] we will discuss how the choice of topology 
on V{E) influences Cb{V{E)). 

The nonlinear dynamics given by the semigroup S^^ is deterministic and 
defines a semigroup of operators on Cb{V{E)), in a way that is reminiscent 
to the equations (I2.ip - (l2.2p but now using the duality structure between 
Ci{V{E)) and V{V{E)). We define, for any G V{E) and ^ G Cb{V{E)), 

rrmn) = ^ (5f^(/in)) . 

The semigroup is called the puUback semigroup of S^^. As we shall see, 
it plays a similar role for the deterministic limiting flow S*^^ associated with 
the equation (II. 2p (lower half of the diagram), as the one the semigroup T/^ 
plays for the flow at the level of the A^-particle system (upper half of 
the diagram): in both cases these are the dual statistical flows. Note that 
for making sense of the pullback semigroup on Cb{V{E)), one needs the 
iiia-P /in ^ S^^{fin) to be continuous, and this will be a major issue when 
all arguments are made precise. 

To connect the A^-particle dynamics with the limiting dynamics, we also 
need mappings between Cb{E^) and Cb{V{E)). On the one hand 

r Cb{V{E)) ^ Cb{E^) 

is dual to the map V{E) 3 ^ ir^f^ £ V{V{E)), and is defined by 

VZGi?^, </)(Z) = (vr^cl.) (Z) = $ (;uf ) , 
where the empirical measure /i^ is defined through (j2.3p . On the other hand 

Cb{E^) ^ Cb{V{E)) 

(j) IH> $ 

is defined in the following way: for each G C},{E^), R^[4>\ is evaluated at 
the point / G V{E) as 

(2.5) f^R^[4>]{f) = R^{f)= [ ^{zi,...,ZN)f{dzi) ■■■ fidzN), 

Jen 

which can be interpreted, as we shall see later, as a real valued polynomial 
taking probability measures as arguments. 



(2.4) 



R 



N 
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2.5. The abstract theorem. We are now ready to give a more precise 
version of the main abstract theorem. The exact statement involves rather 
technical definitions, and its validity depends on five assumptions, (Al) to 
(A5) which are properly stated in Section HI The first assumption is the 
requirement of symmetry under permutations that has already been stated. 
The remaining conditons are (1) estimates on the regularity and stability of 
the nonlinear semigroup, and (2) consistency estimates that quantifies that 
the A^-particle systems and the limiting semigroup are compatible. 

Theorem 2.1 (Fluctuation estimate). Consider a process {Z^)t>o iriE^ /G^, 
and the related semigroups and as defined above. Let /jn G V{E), 
and consider a hierarchy of N -particle solutions = SP{f®^), and a so- 
lution ft = 5*/^^ (/in) to the limit equation. We assume that (Al) to (A5) 
hold. 

Then there is an absolute constant C > and, for any T € (0,oo), there 
are constants Ct,Ct > (depending on T) such that for any N,i G N*, 
with N > 21, and for any 

= (pi (E) ■ ■ ■ (S> ipE jr®^, ^pJ G J^, \\ipj\\^ < 1, 



we have 



(2.6) sup ((5f(/r)-(5f^Un))n,V^' 
[o,r) \^ ^ 

<C^^ + CTf e{N) + Ct i (/in), 

with 

(2.7) n%^{f,^):= [ distg3(^f,/in) /^^(dZ). 

The space T C Cb{E) and the distance distg3 are defined later in Section\^ 
We have used the notation 99 = iS" • • • (8) (fe to denote 
(p{zi,. .. ,ze) = (pi{zi)ip2{z2) ■ ■ ■ iptize), 
and ip 1^-^ to denote 

(99 (g) l^~^){zi, ... ,zn) = (pi{zi)ip2iz2) ■ ■ ■ ipeizg). 

We first note that the left hand side of (|2.6p is the same as the left hand 
side of (|1.3p . the only difference being that the inital data to the A^-particle 
system are assumed to factorize: f^ = f^^, where /in is also the inital 
data to the limiting nonlinear equation. This is a stronger hypothesis than 
merely requiring the initial data to be chaotic (where the initial data to the 
A^-particle system may factorize only in the limit of infinitely many particles) . 
This restriction simplifies the proof, but it can easily be relaxed at the cost 
of some additional error terms. 

The restriction to test functions of the form (pi <Si ■ ■ ■ >pe ]^®(Af-Q^ ^ g_ 
functions depending only on the first i variables, corresponds to analysing 
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^-particle marginals. Hence the theorem imphes the propagation of chaos 
as soon as T is dense in Cb{E) in the topology of uniform convergence on 
compact sets. This condition is satisfied in all examples given below. 

2.6. Main steps of the proof. The proof begins by splitting the quantity 
we want to estimate, 



N-i 



in three parts, each one corresponding to one of the error terms in the right 
hand side of the equation (12. 6p : 



< 



< 
+ 



(2.8) +|(/®^(T^i^^M)o;xf) 



--■■Ti + Ti + Ti. 



Then each of these terms is estimated separately: 

(1) The first term, 7i, is bounded by Cl'^/N^ as proven in Lemma 14.21 
From the definitions of and /x^, it follows that R^[lp\{ijl^) is a 



sum of terms of the form ^p{zj^ 



I , where each index ji , . . . , ji 
describes the set {1, . . . , A'^}. The error is then due to the fraction of 
terms for which two or more of the indices ji, . . . ,ji are the same. 
Hence the estimate of 7i is of purely combinatorial nature, and only 
depends on the symmetry under permutation, i.e. the assumption 
(Al). 

(2) The estimate of the second term, 72, relies on the convergence of the 
A^-particle semigroups to the limiting semigroup Tj°°. This is 
where the A^-particle dynamics and the limiting dynamics are com- 
pared. The estimate can be found in Lemma |4.3[ which depends on 
the consistency assumption (A3) on the generators and the stability 
assumption (A4) on the limiting dynamics. While the generator 

of can be defined in a straightforward manner, defining and es- 
timating the generator G°° of requires a more detailed analysis. 
It indeed involves derivatives of functions acting on V{E), and the 
required differential structure depends on the topology and metric 
structure chosen on V{E). The generator is characterized in 
Lemma l4.ll The assumption (A4) is proved by establishing refined 
stability estimates on the limiting semigroup, showing their differen- 
tiability according to the initial data in a metric compatible with the 
previous steps. 

(3) For the last term Ts, we note that 



R%] (5f^(Mf )) 
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which means that the nonhnear hmiting equation is solved taking a 
sum of Dirac masses as initial data, and an £-fold product of the 
solution is integrated against (p. The resulting function oi Z = 
(zi, . . . , zat) is then integrated against f®^ , which amounts to tak- 
ing an average over all intial data such that the position of the 
particles are independently taken at random from the law When 
A'^ is large, the random empirical measures /i^ are close to /in, i.e. 
np(/if^) ^ in V{V{E)). This implies that the term % vanish 
in the limit when goes to infinity. However the rate of this con- 
vergence sensitively depends on the regularity of the test function 
if and on the continuity properties of S^^ . In all cases considered 
here, the error Ts dominates the other error terms, and effectively 
determines the rate of convergence in the propagation of chaos. The 
precise result, together with the required assumptions, is given in 
Lemma 14.51 

3. Metrics on V{E) and differentiability of functions on 7^(i?) 

This section contains the techical details concerning the space of probabil- 
ity measures on V{E) and its dual, that is the space of continuous functions 
acting on V{E). 

3.1. The metric issue. V{E) is our fundamental "state space", where we 
compare the marginals of the A^-particle density and the chaotic infinite- 
particle dynamics ft through their observables, i.e. the evolution of con- 
tinuous bounded functions on E^ and V{E) respectively under the dual 
dynamics Tf^ and . 

There are two canonical choices of topology on the space of probabilities, 
which determine two different sets Cb(V{E)). 

On the one hand, for a given locally compact and separable metric space 
£, the space M^{£) of finite Borel measures on is a Banach space when 
endowed with the total variation norm: 

yfGM\£), WfWrv ■■= f^iS) + /-{£) 

= sup {f,4>)= sup (/,(/)), 

<l)(iCt(Z), U\\oo<l <l)eCo{£), \\<t>\\oc<l 

where / = — stands for the Hahn decomposition and the equality 
between the two last terms comes from the fact that £ is locally compact 
and separable. 

TV 

We recall that > f (strong topology) when (/fc) and / belongs to 

M^(£) and ||/fe — /||tv ~^ when k ^ oo, and that fk^f (weak topology), 
if 

V99 G Cb{Z) {f,Lp) = lim . 

The associated topology is denoted by a{M^{£), Cb{£)). However, the weak 
convergence can be associated with different, non-equivalent metrics, and 
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the choice of metric plays an important role as soon as one wants to perform 
differential calculus on V{S). 

In the sequel, we will denote by Cb(V{E),w) the space of continuous and 
bounded functions on 'P(-E) endowed with the weak topology, and Ci){V{E), TV) 
the space of continuous and bounded functions on V{E) endowed with the 
total variation norm. It is clear that Cb{V{E),w) C Cb{V{E),TV) since 

TV 

fk > f implies fk^f- 

However, the supremum norm ||$|| loo (75(B)) does noi depend on the choice 
of topology on V{E), and endows the two previous sets with a Banach space 
topology. The transformations tt^ and satisfy: 

(3.1) ||7r^^'||^oo(EiV) < ||^||loo(p(^)) and || J?^[0] ||loo(p(£;)) < ||(/'||Loo(£;iV). 

The transformation vr^ is well defined from Cb{V{E),w) to Cb{E^), but 
it does not map Cb{V{E),TV) into Cb{E^). 

In the other way round, the transformation is well defined from 
Cb{E^) to Cb{V{E),w), and therefore also from Cb{E^) to Cb{V{E),TV): 
for any (j) G Cb{E^) and for any sequence fk so that the weak convergence 
/fe - / holds, we have /f^ - /®^, and then ^ 

The different metric structures associated with the weak topology are not 
seen at the level of Cb{'P{E),w). However any norm (or semi-norm) "more 
regular" than the uniform norm on Cb{V{E)) (in the sense of controlling 
some modulus of continuity or some differential) strongly depends on this 
choice, as is illustrated by the abstract Lipschitz spaces defined below. 

Definition 3.1. Let mg : E — >■ R+ he given. Then we define the following 
weighted subspace of probability measures 

Vg{E):={f eV{E)- (/,m0)<oo}, 

together with a corresponding space of "increments", 

XVg{E) := {/i - ; /i, G VgiE)} . 

If moreover there is a vector space Q (with norm denoted by \\ ■ \\g) which 
contains XVg{E), we then define the following distance on Vg{E) 

V/i, /2 G Vg{E), distg(/i,/2) := ||/i - Mb- 

Remark 3.2. Note carefully that the space of increments IVg{E) is not a 
vector space in general. 

Now we can define a precised notion of equivalence of metrics: 

Definition 3.3. We say that Vg{E) has a bounded diameter if there exits 
Kg > such that 

\^feVg{E), distg{f,g)<Kg 
for some given fixed g &Vg {E) . 
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Two metrics do and di on Vg{E) are said to he Holder uniformly equiva- 
lent on bounded sets if there exists k G (0, oo) and for any a £ (0, oo) there 
exists Ca G (0, oo) such that 

V/i, /2 e BPg,a, ^ [d2{fl,f2)T < dl{fl, f2) < Ca [dsl/l,^)]" 

where 

BPg,a:={ferg{E) ; {f,mg)<a}. 

Finally, we say that two normed spaces Qq and Qi are Holder uniformly 
equivalent (on bounded sets) if this is the case for the corresponding metrics. 

We also define the vector space UC{Vg{E)]'M.) of uniformly continuous 
and bounded function ^ : Vg{E) — )• R, where the continuity is related the 
metric topology on Vg{E) defined by distg above. Observe that this is a 
Banach space when endowed with the supremum norm. 

Example 3.4. With the choice mg := 1, || • \\g := || • \\tv "we obtain 
Vg{E){E) = {V{E),TV) endowed with the total variation norm. 

3.2. Examples of distances on measures when E = R'^. There are 
many ways to define distances on V{E) which are topologically equivalent 
to the weak topology of measures, see for instance [35- 'H]. 

We list below some well-known distances on 'P(R'^) or on its subsets 

P,(M'^) := {/ G P(M^); M,(/) < oo}, q> 0, 

where the moment Mq[f) of order q oi a, probability measure is defined as 

M,{f) := {f , {vY) , {v)^ = l + \v\\ 

These distances are all Holder uniformly equivalent to the weak topology 
a{V{E),Cb{E)) on the bounded subsets 

BPg,a{E) := {/ G VgiM."), Mgif) < a} 

for any a £ (0, oo) and for q large enough. For more informations we refer 

to [12. 



Example 3.5 (Dual-Holder, or Zolotarev's, distances). Denote by dist^; a 
distance on E and fix zq £ E (e.g. zq = when E = in the sequel). 
Denote by LipQ(£') the set of Lipschitz functions on E vanishing at one 
arbitrary point zq £ E endowed with the norm 

mup = mi ■■= sup 

We then define the dual norm: take mg := 1 and endow Vg{E) with 
(3.2) yf,geVg{E), [g-f]l:= sup i^ff^. 
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Example 3.6 (Monge-Kantorovich-Wasserstein distances). For q E [1,(X)), 
define 

Vg{E){E) = V,{E) := {/ G ViE); (/, mg) := (/, dist(., voY) < 00} 
and the Monge-Kantorovich-Wasserstein (MKW) distance Wq by 

(3.3) yf,g€Vq{E), W^{f,g):= inf / distsiz, p{dz,dz), 

where Il{f,g) denote the set of probability measures p G 'P{E x E) with 
marginals f and g (p{A x E) = f{A), p{E x A) = g{A) for any Borel set 
A C E). Note that for Z, Z ^ E^ and any q G [1, 00), one has 
(3.4) 

Wg (^fJ-Z^fJ-^) = di,(EN/6M)i^^^) ■= ( ^5^disti7;(zi,Z<^(i))'?J , 

and that 

(3.5) yf,gePiiE), Wi{f,g) = [f-g]l= sup {f-g,4>) 



as well as 



(3.6) yq 00), yf,gerg{R''), Wi{f,g)<Wg{f,g). 

We refer to [41\ and the references therein for more details on the Monge- 
Kantorovich- Wasserstein distances and for a proof of these claims. 

Example 3.7 (Fourier-based norms). For E = W^, mg := 1, let 

yfeTVg{E), ||/||g = |/|, := sup^, s > 0. 

We denote by (which includes XVg{E) for s large enough) the Banach 
space associated to the norm \ ■ \s- Such norms first appeared in connection 
with kinetic theory in |16| . 



Example 3.8 (Negative Sobolev norms). For E = W^, mg := 1, let 



yfeTVgiE), ||/||e = ||/||H- 



, s>0. 



We denote by (which includes XVg{E) for s large enough) the Hilbert 

space associated to the norm \\ ■ 

For E = M'^, mg := 1, and some integers k,i > 0, we also define 

V/ G TVg{E), \\f\\g = \\f\\H-UR'^^ ■= ^up if,^), 



where 
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We denote by H 



^) (which includes ZVg{E) for k large enough) the 



Hilbert space associated to the norm 



3.3. Differential calculus for functions of probability measures. We 

start with a definition of Lipschitz regularity, for which a mere metric struc- 
ture is sufficient. 

Definition 3.9. For metric spaces Qi and Q2 we denote by C^'^{Qi,Q2) the 
space of functions from Qi to Q2 with Lipschitz regularity, i. e. the set of 
functions ^ : Qi ^ Q2 such that there exists a constant C > so that 

(3.7) V/,5Gai, dist^^(^'(<7),^(/)) <Cdist^^(g,/). 

We then define the semi-norm [•](^o,i(g^ g^-^ on C^''^{Qi,Q2) cis the infimum of 
the constants C > such that (13.71) holds. 



Tlie next step consists in defining a higher order differential calculus; this 
is where the assumption that metrics are inherited from a normed vector 
space structure plays a role. 

Definition 3.10. Let Qi and Q2 be normed spaces, and let Qi and Q2 be two 

metric spaces such that Gi — Gi C Gi- For A; G N, we define C^'^{Gi', G2) to be 
the set of bounded continuous functions ^ : Gi ^ G2 such that there exists 
: Gi — B^{Gi,G2) continuous and bounded, where B^{Gi,G2) is the 
space of bounded j -multilinear applications from Gi to G2 (endowed with its 
canonical norm) for j = 1, . . . ,k, and some constants Cj > 0, j = 0, . . . , k, 
so that for any j = 0, . . . ,k 



(3.8) \ff,geGi, 



^{g)- Y.{D^^{f),{g-fY 



1=0 



<C,\\g-f\\i: 



Q2 



(with the convention D^"^ = 

We also define the following seminorms on C^'^{Gi,G2) 

[%o:= sup ||I)%(/)||„,,. .7 = 1, 

/eai 



with 



and 



\L\ 



»(Gi,G2) 



sup 



sup \\L {hi, . . . , hj) 

h^,\\h,\\g^<l,l<i<j 



n9)-EU{D'^if)d9-fr 



k, 



102 ' 



02 



Finally we combine these semi-norms into the norm 



k 

E 



k,l- 
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Remark 3.11. Observe that for any j > 1 the LHS of (13. Sp makes sense 
since 



1=0 



=1 

since ^-(5) - *(/) £ G2-G2 C G2 and D^^{f) G ^^'(^1,^2)- Then note that 
our definition is very close to the usual Frechet definition of differentiability 
in Banach spaces for the function h ^{f + h) with h = g — f £ Gi, except 
that the domain and range are restricted to subsets that have no vectorial 
structures and are not open within Qi and Q2- We also only consider Lipschitz 
differentiability. 

The following lemma confirms that this differential calculus is well-behaved 
for composition, which seems to be a minimal requirement for further appli- 
cations. 

Lemma 3.12. Consider U G C^'^{Gi,G2) and V G C'='H^2,^3)- Then the 
composition := VolA belongs to C^'^{Gi,Gz)- Moreover the following chain 
rule holds at first order k = 1 

(3.9) V/G^i, D^[f] = DV[U{f)\oDU[f], 

with the estimates 

' [*]o,i < [V]o,i [Z^]o,i, 

[*]i,o < [V] 1,0 [Z^] 1,0, 

, [*]i,i<[V]i,o[^]i,i + [V]i,i[W]2,. 

At second order k = 2 one also has the chain rule 
(3.10) 

V / G ^1, DH[f] = D^Uif)] o {DU[f]®DU[f]) + DV[U{f)]oD^U[f]. 

Proof of Lemma l3. J2[ It is straightforward by writing and compounding 
the expansions of lA and V provided by Definition 13.101 □ 

3.4. The subalgebra of polynomials in Cb{V{E)). In Section [2] we de- 
fined a map : C{E^) — )■ Ch{V{E)), which may be used to define a subal- 
gebra of polyomials in Cb{V{E)). We first define the monomials: 

Definition 3.13. A monomial in Ci,{V{E)) of degree i is a function R^[(p] 
with ip = ipi ■ ■ ■ (p£, ipi G Ch{E) and £ G N. Explicitly 

R%]{f) = [ ip{zu...,ze)df^\zi,...,ze) 

e 



n / ^i{z)df{z), 

i=i 
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which is well-defined for all f G 'P{E). 

The product of two monomials is defined in a natural way 
^^i+^2j^0^j^ and the polynomial functions are linear combinations of mono- 
mials. These form a subalgebra of Ch{V{E)) that contains the constants and 
separates points in 'P{E), and hence the Stone- Wcicrstrass Theorem implies 
that this subalgebra is dense in Cb{V{E)), where the meaning of "dense" 
depends on the topology chosen on 'P(-E). 

While the polynomials in R always are differentiable, the smoothness of the 
polynomials depends on the metric structure. We need first some preliminary 
definitions. 

Definition 3.14. • Duality of type 1: We say that a pair {T,Q) of 
normed vector spaces such that J- C Ch{E) and V{E) — V{E) C Q 
satisfies a duality inequality if 

(3.11) S^f,geViE), V(^G \{{f-g),^)\<C\\f-g\\gMr. 

• Duality of type 2: More generally we say that a pair {J^,Vg{E)) of 
a normed vector space T C C\j{E) endowed with the norm \\ • \\jr and 
a probability space Vg{E) C 'P{E) endowed with a metric dg satisfies 
a duality inequality if 

(3.12) ^ f,gerg{E),^^e:F, \{g - f , ^)\ < C disig{f , g) Mr . 

Lemma 3.15. If ^ and the pair {J^,G) satisfy a duality of type 1, the 
polynomial function R^[(p\ is of class C^'^{Vg{E),'E.) for any k > 0. In the 
more general case where the pair {J^,Vg{E)) satisfies a duality of type 2, the 
polynomial function R^l^p] is at least of class C^'^{Vg{E),'K). 

Proof. It is clearly enough to prove the lemma for monomials, and the proof 
then mainly follows from the multilinearity of R. In the case of duality of 
type 2, then the conclusion follows from 



R%]{f2)-R'[mi) = E n (^«,/2-/i) n^^^'/i) 

1=1 \l<fe<i / \i<k<e 

In the case of duality of type 1, we define 



g^K, DR'[^]{f){h) ■.= ^ll[{^j,f)] {^^,h), 

1=1 
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and we write 
R'M{f2) - R%]ifi) - DR%]{h){h - h) = 

= E ( n ^^kj2)\ {^jj2-fi) I n i^kji)] X 

l<j<i<e \l<k<j J \j<k<i J 

X {'^i,f2 - /l) n ("PkJl) 



,i<k<i 



We deduce then 

R%]{f2) - R'Mfi) < \Mi,T^(L^r-^ 11/2 - /lib, 



DR%]ih)ih) 



< 



R%]{f2) - R'Mfi) - DR%]{h){h - h) < llv'lli,^2^(Loo).-2 11/2 - fi 



where 

ll¥'lll,J-fc®(L°°)^-fc • = 
< 

and we have defined 

\\v\\oo,T''»(L°°y-'^ • = 

< 



il,...,ik distinct in [\l,i\] iy^Cu v ■ -i**;) 



l¥'lloo,J-2^{L°°)<-2 



for A; = 2, 



max llv'jillj" • • • llV'iJb IT 

ll,...,lk distinct in [| 1,^1] . , . , 

J^(il,...,Jfe) 

llv^l' 



L°°{E) 



since || • \\lo-^e) < II • lb- This proves that i?^^ G C^'^{rg{E),M.). The 
cases k > 2 are proved similarly. □ 



4. Assumptions and technical lemmas 

In this section we collect the lemmas used in the proof of Theorem 12.11 
and the technical assumptions that are needed. 

The first assumption is simply the statement that the A^-particle dynamics 
is well defined and invariant under permutation of the particles. 



(Al) A^-particle semigroup. The family consists of strongly 
continuous Markov semigroups on Cf){E^) that are invariant 
under permutations of indices. We denote their generator 
by and we denote by the dual semigroup on the A^- 
particle distributions. 
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4.1. The generator of the pullback semigroup. While the definition of 
and its generator is rather standard, it takes more care when defining 
the puhback of the nonhnear semigroup and the corresponding generator. 
The second assumption that we need to impose on the system is related to 
this, and relates to our definition of a differential calculus of functions on 
V{E). 



(A2) Nonhnear semigroup. Consider a probability space 
Vg^ (E) (defined in Definition 13. ip associated to a weight func- 
tion mgj, endowed with the metric induced from a normed 
space Qi, and with bounded diameter. Assume that for any 
r > we have: 

(i) The equation (jl.2p generates a continuous semigroup 
S[^^ on Vg-^{E) which is uniformly Lipschitz continuous: 
there exists Cr > such that 

V /, 5 e rg, (E), sup distg, < a distg, (/, g). 

tG[0,T] 

(ii) There exists 6 £ (0, 1] such that the (possibly nonlinear) 
generator Q (introduced in equation (|1.2p ) is bounded 
and 5-Holder continuous from Vg^{E) into Qi in the 
following sense: there exist L,K > so that for any 
f,geVg,{E) 

\\Qif)\\g,<K, \\Q(f)-Q(g)\\g^<L\\f-g\\'g^. 



This assumption is sufficient for defining the generator of T^: 

Lemma 4.1. Under assumption (A2) the pullback semigroup is a con- 
traction semigroup on the Banach space [/C(Pg^ (i?); M) and its generator 
G°° is an unbounded linear operator on UC{Pg^ {E);M) with domain Dom(G°° 
containing C^'^{Pg-^^{E);W). It is defined by 
(4.1) 

e C'^HPg,{E);R), V/ G Pg.iE), (G^^) (/) = (Z?<I>[/], Q(/)) . 
Proof. The proof is split into several steps. 

Step 1. We claim that for any /;„ G Vg^ (E) and r > the map 

S{f,,,):[0,T)^Vg,{E), i^5f^(/,„) 

is right-differentiable at t = 0+ with 5(/in)'(0) = Q{fin)- 

Denote ft := S'f^/in- First, since Q{ft) is bounded in Qi uniformly on 
t G [0,r] from (A2)-(ii), we have, uniformly on /in G Vg-^{E)^ 



(4.2) 



Wft — /inllei 



Q{fs)ds 



< Kt, 
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and then using (A2)-(ii) and the inequahty (j4.2p we obtain 

Wft-fin-tQiUWg, = [\Q{fs)-Q{fin))ds 

Jo 



< 



L / ll/s - /inllg, ds 

Jo 

ft f 
L / {Ks)^ds = LK^- 
Jo 1 



5' 



which imphes the claim. 

Step 2. We claim that {T^) is a Co-semigroup of linear and bounded (in 
fact contraction) operators on UC{Vg-^^{E);M.). Indeed, first for any <1> E 
U C {Vg^{E)]'K) and denoting by uo^ the modulus of continuity of we have 

\{Trm)-{Trm)\ = \Hs^\g)) - ^{s^\9))\ 

< a;*(diste,(5f^(5),5f^(/))) 

< LJ.s>iCr distg^{ f,g)) 

so that T^^ E UC{Vg^{E);R) for any t > 0. Next, we have 

\\T^°°\\ = sup llTj^^ll = sup sup |$(5f^(/))| < 1, ||$|| = sup |$ 

II*II<1 II*II<1 /GPSiC-B) h£Vg^(E) 

Finally, from (gj]), for any $ E [/(^(Pg;, (-E); M), we have 

\\T^^-^\= sup |$(5f^(/)) - «>(/)| < ^ ast^O+. 

/eT'g, (E) 

As a consequence, Hille-Yosida Theorem (see for instance |331 Theorem 3.1]) 
implies that (Tj°°) is associated to a closed generator G°° with dense domain 
dom(G~) C UC{Vg,{E)-m). 

Step 3. A candidate for this generator is defined as follows. Let G°° be 
defined by 

E C''HVg,{E);R), V/ E Vg,{E), (G~cl>)(/) := {D^fiQif)) ■ 

The RHS is weh defined since D^{f) E B{gi,R) = g[ and Q(/) E ^i by 
assumption. Moreover, since both / i— D^[f] and / i— t- Q{f) are uniformly 
continuous so is the map / ^ (G'°°$)(/). It yields (5°°$ E UC{rg,{E);M.). 
Step 4. Finally, by composition, 

yfGVg,{E){E), t^TrHf) = ^oS^'^{f) 

is riglit-differentiable at t = 0^ and 

^(rri>)(/)k=o := ^($o5(/)(t))|t=o 

I)ci>(5(/)(0)),^5(/)(0) 



(D<I.[/],Q(/))= G°°<I> (/) 
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which impHes $ G Dom(G°°) and (14. ip . 



□ 



4.2. Estimates in the proof of Theorem 12.11 The proof of Theorem 12. II 
rehes on the three following lemmas, together with their assumptions. 



4.2.1. Estimate of 71. The first error term is estimated by the following 
combinatorics argument. 

Lemma 4.2. Let , and R^[^] as defined in Section\^ (see FigureU^), 

and let ip £ C^{E'^). For any t > and any N >2l 

(4.3) 



Ti 



Proof. Since S^{f^^) is a symmetric probability measure, estimate (14. 3p is 
a direct consequence of the following estimate: For any if G Cb{E^) and any 
iV > 2£ we have 



(4.4) 



< 



N 



Here the symmetrized version of a function (j) G Ci,{E^), is defined as 

(4-5) 4ym = Tg^ (f>a- 



N\ 



o-G6jv 



As a consequence for any symmetric measure / G Vsym{E ) we have 



(4.6) 



{f'',R%]{^^^^))-{f^\ip) 



pN 



< 



N 



To establish the inequality (j4.4p . we let i < N/2 and introduce 
AN,e:={{ii,...,ii)£{l,...,NY : yk^k', ik^ik'} and BN,e ■= A' 

Since there are A^(A^ — 1) . . . (A^ — ^ + 1) ways of choosing £ distinct indices 
among {1, . . . , A^} we get 



Ni- 



ls 



1- 1 



1 

N 



... 1 



1 



£-1 



< 1-exp -2j^ 



A^ 



1 — exp In I 1 



N 



i=0 



N j - N' 



where we have used 

VxG [0,1/2], ln(l-x)>-2x and Vx G M, e""^ > 1 - x. 
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Then we compute 

1 ^ 

ji,...,i^=l 

{h,...,ie)€AN,e (ii,...,jf)GBjv/ 
1 1 . . /i' 



(N -£)l 

1 v-^ /2£'^ \ 



o-e6jv 

and the proof of (j4.4p is complete. Next for any E V{E^) we have 

\ \ / sym / 

and (gH) follows from (g^]). □ 

4.2.2. Estimate o/Ti- The second error term is estimated thanks to a con- 
sistency result for the generators of the particle system and the limiting 
dynamics, and stability estimates on the limiting dynamics. To proceed we 
need to introduce the two corresponding assumptions. 



(A3) Convergence of the generators. Let Vg^ [E) be the probabil- 
ity space endowed with a metric considered in (A2). For k = 1 
or 2 there is a function £{N) with lim7v-s>oo e(-^) = such that 
for ah iV e N: and ah $ G C'^'^Pg, (^); M), the generators 
satisfy 

(4.7) ||G^(7r^ci,)_vr;vG~(cI>)||^^(^,) <e(iV)||$||^M(p,^(ij)). 

where Q is the nonlinear operator involved in equation (II. 2p . and 
G°°(<J>) = {Q,D^) is the generator of the pullback semigroup 



(A4) Differential stability of the limiting semigroup. We assume 
that the flow S^^ is C'^'^iVg^ {E),Vg^ (E)) in the sense that there 
exists Ct > such that 

(4-8) / ||'S'f^||r'fe■l^p^.^B^p^_^R^^ < Ct 



\Ck^[Vg^(E),Vg^[E)) 

where Vg^ (E) is the same subset of probabilities as Vg^ (E) , en- 
dowed with the norm associated to a normed space Q2 D Gi 
possibly larger than Gi and where k is the same integer as above. 
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Lemma 4.3. Suppose that the assumptions (Al) to (A4) are satisfied, and 
let T2 be as defined in equation Ii2.8\) . Then 



(4.9) T2 := |(/r,7^*^(«'M°/^f))-(/r,((rr^'M)°/^f 

for an explicitly given constant C{k,£) depending only on k and i. 
Proof. We start from the following identity 

Jo 

which we evaluate on <I> = R^[cl)\ G ChiJ-'iE)). From assumption (A3) we 
have for any t G [0,7"] 



Tl'TlNR'W]-TTNTrR'M 
t 



S!L, (/^) , [G'^TT^ - TT^G^] {TrR%])) ds 



< 



[ [G^'ttn-ttnG^] {T^R'[ip]) 
Jo 



ds 



(4.10) 



T^R%] 



ds. 



Since T^^^ (i?^ [(^] ) = R^[tf] o with E C'''\Vg,{E),rg,{E)) thanks 
to assumption (A4) and R'^[ip] G G^'^{Vg2{E),M.) because tp G jf^ (see 
Subsection 13. 4p . we obtain with the help of Lemma |3. 121 that Tj°° (ii^ [99] ) G 
G^'^{Vg-^{E)) with uniform bound. We hence conclude that 



(4.11) 



ds<C{k,e)GT R^[ip] 



where C{k, I) < l"^ since A; = 1 or A; = 2. 

Going back to the computation (|4.10p . and plugging (|4.1ip we deduce 

m. □ 



4.2.3. Estimate of T3. The error term T3 from equation (12. Sp depends on 
estimating how well the initial data for the nonlinear equation /in can be 
approximated by an empirical measure, and how well this error is then prop- 
agated along the semigroup. To this purpose we need to make a stability 
assumption for the limiting semigroup. 
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(A5) Weak stability of the limiting semigroup. There is proba- 
bility space Pg^{E) (corresponding to to a weight function mg.^ 
and metric distgg) such that every T > there exists a constant 
Ct > such that 

(4.12) 

V/i,/2 G PgsiE), supdistg3 (5f^(/i),5f^(/2)) < Ct distg.ih, h). 

[0,T) 



Remark 4.4. Observe that when Vg^{E) = Vg^{E) = 'Pg^{E) with the same 
weights and distances, the assumption (A5) is included in (A4). However 
it is crucial to have the flexibility to play with different metric structures in 
these assumptions. 

Lemma 4.5. Assume that the limiting semigroup Sf'^ satisfies assumption 
(A5) for some probability space 'Pg^{E). Let satisfies a duality inequality 
with Vg,^{E) as defined in Definition \3.14\ Let ip G Ci){E^). 
Then for any G V{E), t > and N > 21 we have 

(4.13) Ts := |(/r, (rrii^M) o^f ) - ((Sr(/J)^\^) 

< ICT^^ifin) ll¥'llj-3(g.{L°°)^-l, 

where Q,^{fia) is defined in (|2.7p and \\'p\\j^^^(L°°y-'^ defined in the proof 
of Lemma \3.15[ 

Proof. We split Ts in two terms, the first one being 

7-3,1 := (fi^{TrR%])o^,f 

R%] {S^il^^)) fuM^i) ... fUdzN) 



j^^ {^^a^{Z)^ /in(dZl) ... /in(dZ7v), 



with 

Vi = 1, . . . ai = a,{Z) := f p,{w) 5f^(^f )(du;). 

J E 

Similarly, we write for the second term 

r3,2 = ((5f^(/in))''',V')=y^^(^n6^j /i„(dZl).../i„(dz;v), 

with 



Vi = 1, . . . h:= ^,{w) S'^'^{f,^){dw) 



E 



Using the identity 

til 

Oi - 6i = ^ oi . . . Oi-i (oj - hi) hi+i ... 



1=1 i=l 1=1 
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we get 



(4.14) Ts<Y,\IIM\l^^ 



i=l 




ai{Z) - hi\ /in(d2i) . . . fi^{(izN)- 



Then by using the duahty bracket together with assumption (A5) we have 



\ai{Z) - bi 



< 11^,1^3 distg3(5f^(A0,5r(/xf)) 

(4.15) < C'T||(/?i|l^3distg3 (/i„,/xf). 

Therefore combining (14.14^ and (14.15^ (for any 1 < i < ^), we conclude that 
KWi holds. □ 

In order to use Lemma [4.5l we need an estimate on the term (/in). This 
information is provided by the following quantitative version of the law of 
large number for empirical measures taken from |35] . We refer to |32j for a 
more detailed discussion of this issue. 

Lemma 4.6. For any /;„ G Vd+5(^'^) and any N > 2 there exists a constant 
C which only depends on d and M^+^^f in) so that 



(/in) 



W^2(;uf,/,0Vr(d^)<CiV-^. 



\2 f^N 



4.3. A remark on assumption (A4). In this section we briefly explain 
how our key estimate (A4) can be obtained in the case of a nonlinear oper- 
ator Q which splits into a linear part and a bilinear part: 

(4.16) V/GP(S), Q{f) = Qiif) + Q2{f,f) 

with Qi linear and Q2 bilinear symmetric. 

For two initial data /in and g^^ in a space Pg{E) of probability measures, 
and some initial data /lin G we introduce the following evolution equations, 

{ dtg = Q{g) = Qi{g) + Q2{g,g), g\t=o = 9in, 

dtf = Qif) = QlU) + Q2(/, /), f\t=0 = fin, 

[ dth = DQ[f]{h) = Qi{h) + 2Q2{f,h), h\t=^ = hn 

In the third equation the solution ht depends linearly on h^^ (but also 
nonlinearly on /in): it is formally the first-order variation of the semigroup, 
i.e. DS^'^fiM = h. 

We now want to write a second-order variation of the semigroup. To this 
purpose we conside ht and ht two solutions to the third equation above and 
write 



dtr = DQ[f]ir) + -D-'Q{f){h,h) = Qi{r)+2Q2{f,r)+Q2{h,h), r^t=o = 0- 
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In this fourth equation rt depends bihnearly on ht and ht which are two 
solutions to the third equation, and therefore bihnearly on /iin,/iin (it also 
depends nonlinearly on again): it is formally the second-order derivative 
of the semigroup D'^S^^[fin]{hin, /ii„) = r^. Observe that the initial data Tin 
are always zero for this second variation problem since the map {ft)t>o /in 
is linear. 

Consider now hi^ = /lin = ^in — /in (which implies ht = ht). Let us define 
s := f + g, d := g — f , U! := g — f — h, := g — f — h — r, for which we get 
the following evolution equations 

' dtd = Ql{d) + Q2(s,d), d\t=0 = din = g-m - /in, 

< dtu = Qi{uj) + Q2{5,uj) + Q2{h,d), uj\t=o = 0, 

. dtip = QiW + Q2{s,i;) + Q2{h,uj) + Q2{r,d), V'|t=o = 0. 

Now we can translate the regularity estimates on Sf'^ in terms of estimates 
on these solutions on some given time interval [0, T]: 

f sup lldtlb, < Crildinlle, =^ Sf"^ €C''HPg,iE),Pg,{E)) 

iG[0,T] 

' [ SUPtg[o,T] Iktllea < CTlldinll^^ J 

r S^^eC'^\Pg,{E),Pg,{E)) j 

< sup,e[o,T] llnllg. < CT||/iin||gJ|/iin||e, > ^ G C''^(Pe,(i?), Pg,(i?)). 

. ( SUPtg[o,T] Ut\\g2 < C'Tlldinllg^ J 

Such estimates are typically obtained by energy estimates for the equations 
satisfied by d, r, cj and ip, for a well chosen "cascade" of norms connecting 
II ■ Ibi to II ■ llfo (s6e later in the applications). 



5. Maxwell molecule collisions with cut-off 

5.1. The model. In this section we assume that E = M'^, d > 2, and 
we consider an A^-particle system undergoing a space homogeneous random 
Boltzmann collisions according to a collision kernel b G L^([— 1,1]) only 
depending on the deviation angle and locally integrable. This is usually 
called Maxwellian molecules with Grad's angular cut-off, as introduced in 
[23| [2^ [29] . We make the normalization hypothesis 

11611,^1 = / b{a^)da=l. 

Let us now describe the stochastic process. Since the phase space i?^ 
corresponds to the velocities of the particles, we shall denote Z = V in 
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this section. Given a pre-collisional A'^-system of velocity particles V = 
{vi, . . . ,vn) € = (W^)^ , the stochastic runs as follows: 

(i) for any i' ^ j' , we draw randomly for the pair of particles (v^/jVjr) 
a random time Ti'j' of collision according to an exponential law of 
parameter 1, and then choose the collision time Ti and the colliding 
pair {vi,Vj) (which is a.s. well-defined) in such a way that 

Ti = Ti,j := min Tj/ 

(ii) we then choose a G at random according to the law b{cos9ij) 
where we define the angular deviation O^j by cos 6ij = cr-(vj —Vi)/\vj — 

(iii) the new state after collision at time Ti becomes 

V* = V*j = Rij,aV = {vi,...,v*,....,v*,...,vn), 
where the rotation Rij^a on the {i, j) pair with vector a is defined by 

with 

Wij = Vi + Vj, u*j = \uij\a, Uij = Vi- Vj. 

Scaling the time by a factor 1/N and repeating the above construction 

lead to the definition of a Markov process (Vf ) on (W^)^. It is associated to 
a Feller semigroup (T^^) with generator . Moreover the master equation 
on the law is given in dual form by 

(5.2) ^(/*^^) = (/f,G^<^) 
with 

1 ^ r 

(5.3) (G^^)(V) = - Yl / b{cos9^j)[^*j-^]da 

l<i<j<N-^^'' ^ 

where ip*j = ^{V*j) and if = ^{V) G ^(M^'^). Finally, the fiow f-^ ^ //^ 
defines a semigroup for the A?^-particle distributions which is nothing but 
the dual semigroup of T/^. 

Note that the collision process is invariant under permutation of the veloc- 
ities, and satisfies the microscopic conservations of momentum and energy 
at any collision time 

TV 

Va = l,...,d, Y.''ka = J2^ka, |^*|' = |V^|':=El^*^l'- 

k k k=l 

We write V = {vi)i<i<N = (vi, ■ ■ .,vn) G E^^ and v = {va)i<a<d G IK'^, 
so that V = {via) e R^'^ with Via G M. 
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As a consequence, for any symmetric initial law £ VgymiM^'^) the 

law density remains a symmetric probability and conserves momentum 
and energy 

Va = l,...,d, / (yvka^fl'idv)-- 



V0:M+^M+, / e{\V\')ff{dv)= e{\v\')fi''idv 





The formal limit of this A^-particle system is the nonlinear homogeneous 
Boltzmann equation on V{W^) defined by 

(5.4) §-^ft = Qift,ft) 

where the quadratic Boltzmann collision operator Q is defined by 

(5.5) {Qif, f),ip) := [ bie) {cl>{wl) - <I>{W2)) da f{dwi) fidw2) 
for ip £ CbiR"^) and / G V(R'^), with 

* WI+W2 , \W2-Wi\ ^ WI+W2 |u^2-'U^l| 

(5.6) w^ = + a, W2 = a 

and cos^ = a ■ {v — w)/\v — w\. The equation (I5.4p - (l5.5p is the space homo- 
geneous Boltzmann equation for elastic collisions associated to the Maxwell 
molecules cross section with Grad's cutoff. We refer to the textbooks [lOj and 
|40) and the numerous references therein for both the physical background 
and the mathematical theory of the Boltzmann equation. This equation gen- 
erates a nonlinear semigroup S^^ on 'P(M'^) defined by S^^ fin '■= ft for any 
/in G T'{W^), which satisfies conservation of momentum and energy: 

yt>0, [ vft{dv)= [ vfUdv), [ \v\^ft{dv)= [ \v\^fUdv). 

JW jRd 

5.2. Statement of the result. On the one hand, it is well known that for 
the collision kernel that we have chosen the A^-particle Markov process (V/^) 
described above is well defined for any initial velocity , and in particular, 
for any given initial law S 7^syin((]K'^)^) there exists a unique solution 
jN g 'Pgyjjj((M'^)^) to equations (|5.2p - (j5.3p so that the A^'-particle semigroup 
is well defined, see |24| |251 l38l l30j . On the other hand, it is also well 
known that for any /;„ £ 7^g(]R'^), q > the nonlinear Boltzmann equa- 
tion (I5.4p - (l5.5p has a unique solution ft G Vq{M.'^). This solution conserves 
momentum and energy as soon as q >2, see for instance [38| [39l [T^ I40| . 

Our mean field limit result then states as follows. 

Theorem 5.1 (The Boltzmann equation for Maxwell molecules with Grad's 
cut-off). Consider an initial distribution f^ G 'Pq{M.'^), q > 2, the hierarchy 



30 



S. MISCHLER, C. MOUHOT, AND B. WENNBERG 



of N -particle distributions //^ = S[^{f^^) following (|5.2p . and the solution 
ft = 3^"^ if,,,) following 

Then there is a contant C > and, for any T > 0, there are constants 
Ct, Ct > such that for any 

<^ = ipi (g) ■ ■ ■ (g) ifi £ J" := C;,(R'^) nLip(R'^), H^JjUj- < 1, 

we have for N > 2i: 
(5.7) 



sup 

[0,T] 



where $7]^^ was defined in (|2.7|) and W2 is the quadratic MKW distance 
defined in (13. 3|) . 

As a consequence of (|5.7|) and Lemma \4-(>\ this implies propagation of 
chaos with rate e{N) < C{£, T, N^^+i for any initial data /i„ G Pd+5(M'^), 
where C{£,T, f-m) is an explicitly computable constant. 

For the Boltzmann equation with bounded kernel, propagation of chaos 
has been estabhshed by McKean in [29], where he adapted the method intro- 
duced by Kac in [24] based on the Wild sum representation of the solutions to 
the Boltzmann equation. Griinbaum in [21j gave an alternative proof based 
on the same "duality viewpoint" as developed in the present paper. Sznit- 
man in [36] also gave a proof of propagation of chaos based on a nonlinear 
martingale approach. In all these works, propagation of chaos is proved but 
without any rate of convergence (as the number of particles goes to infinity). 
Graham and Meleard in \18\ [T9| [30] were then able to prove the propagation 
of the chaos with the sharp rate C{i, T) /N . Their proof is based on the con- 
struction of a stochastic tree associated to the process which is specific 
to the Boltzmann equation with bounded kernel. More recently Kolokoltsov 
in [26] proved a fluctuation estimate for similar processes using a "duality 
view point" like the one developed by Griinbaum and used also in our work. 
His fluctuation estimate is similar to our Lemma [4.9l (and of the same order), 
but pays less attention to the remaing terms of our estimate. Fournier and 
Godinho [ 15] prove the propagation of chaos for a one-dimensional caricature 
of the Boltzmann equation using a coupling method in the spirit of [381 137] . 
Their chaoticity estimate is of the same rate as ours. 



5.3. Proof of Theorem [STU The assumptions (A1)-(A2)-(A3)-(A4)- 
(A5) needed to apply Theorem 12.11 will be verified step by step. In this proof 
we fix = = Cq{W^) and J"3 = Lip(M'^) and define Vg^{E) = Vg^{E) := 
ViW'') endowed with the total variation norm || • \\tv, 'Pg^iE) ■= P2(M'^) 
endowed with the quadratic MKW distance W2- Notice that {Qi,J-i) and 
(^2) -^2) satisfy a duality inequality of type 1, and {Pg^,J''i) satisfy a duality 
of type 2 (see Definition I3.13P . 
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Proof of (Al). The symmetry assumption is satisfied because of the weU- 
known properties of the Bohzmann-Kac A^-particle system, and we refer to 
the previous works [29l EH [l8l [H [30] for details. 

Proof of (A3). We claim that there exists Ci G E+ such that for all 
(5.8) 

which is nothing but (A3) with k = r] = 1 and £{N) = Ci N~^. 
Take $ G C^'\Vg^{E),R), set = D<^[fi^] and compute 



- y 



l<i<j<N 
l<i<j<N 

On the one hand, we have 

N 



da 



O ( ||^>|bi.i 



11^ - 



TV 



i=hiV)) 



da {=I2{V)). 



'1 



-y 



2A^2 

1 
2 



6(%) [</.«) + cPiv*) - cl>{vi) - cl>{vj)] da 



b{9) mv*) + (^{w*) - <P{v) - ^{w)] fi^idv) {dw) da 



= (Q(/x^,/i{y),<A). 

On the other hand, we have 



h{V) 



< 




Collecting these two terms we have proved that (jS.Sp holds. 

Proof of (A4). Here we prove that for any /, /i G P(M'^) and for any 
T > 



(5.9) sup S^'^{g)-S^\f)-CST[n{9-f) 

*G[0,T] 



< e4||7|| 



TV 



\9-f\\ 



where CS'^[f] is the linearization of 5^ at /. As a consequence, this implies 
that ( A4) holds with k = r] = 1 and the previous definitions of Vg^ {E) and 
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Vg2{E). We denote by ft,gt, ht the solutions to the following equations: 

dtft = Qift, ft), f\t=0 = fin, 

dtgt = Q{gt,gt), g\t=o = ffm, 

dtht = 2Q{ft, ht) := Q{ft, ht) + Q{ht, ft), h\t=o = /ii„, 

where Q denotes the symmetrized form of the bilinear collision operator. 
The third equation corresponds to the first-order variation of the semigroup: 
£5'^[/](/iin) = ht solution to this equation. 

Standard Gronwall arguments show the existence and uniqueness of such 
solutions, which moreover satisfy, uniformly on [0, T] 

2 T 2 T 

WhtWrv < e WhinllTV, \\gt — ftWrv < e llftn — /in||TV- 

Next, writing rt := gt — ft — ht, we find that this expression satisfies the 
equation 

dtn = Qift + gt,rt) + Q{gt - ft, h), = 0. 

Introducing yt := ||?^t||Ty, we have 

y't < l\\Q{ft + gt,rt)\\Tv + \\Q{gt-ft,ht)\\Tv 

< II7II00 Wft + gtWrv WnWrv +C\\gt - ftWrv WhWrv 

< Cyt + Ce^'\\h-f\\ly, 

from which we deduce 

VtG[0,r], yt < e''^ \\h - ffTV ■ 
This concludes the proof of (I5.9P . 

Proof of (A2). Assumption (A2)-(i) is clearly a consequence of (A4). For 
(A2)-(ii) we write 



\\Qif,f)-Q{g,g)\\Tv = sup [ {Q{f,f) - Q{g,g))^dv 

\M\l°^ <i Je 

= sup / {ff*-gg*) / b {ip' - If) da dv dv:, 

\\ifi\\LOc<l J ExE Jgd-i 
< m\L^\\f-9\\TV, 
so that the function / Qif,f) is Lipshitz from Vg^{E) to M'^(M'^). 

Proof of (A5). It is known since the seminal work of Tanaka |38| that 
the nonlinear Boltzmann fiow associated to Maxwellian molecules is a con- 
traction for the quadratic MLW distance W2- for any g^a G Vi{M.'^) the 
solutions ft, gt to the Boltzmann equation (j5.4p satisfy 

sup W2{ft,fgt)< W2{f in, g-n.)- 

[o,r] 

That immediately implies (A5) in the space Pg^{E) defined above. □ 
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6. VLASOV AND McKeAN-VLASOV EQUATIONS 

6.1. The model. In this section we assume that E = (where m = d or 
m = 2d with d the physical space dimension, see later) and we consider an A^- 
particle system which undergoes McKean-Vlasov type stochastic dynamics, 
i.e. a drift deterministic force field combined with diffusion. We refer to 
the lecture notes \37\ |30] and the references therein for more details on the 
model, and among many references, we highlight the recent paper [5] for 
recent results and references (using the so-called "coupling" method). The 
method we shall present here does not rely on any of these references. The 
results in this section are mostly not new but compare to the latest results 
of mean field limit on this equation as far as we know. Indeed we shall make 
strong smoothness assumptions on the coefficients of the evolution equation 
in order to avoid technical difficulties and our goal is to advocate for our new 
method and show its power and ability to deal with very different models. 

We assume that the A'^ particles = (-Zi,t, . . . ,Zjq^t) satisfies the sto- 
chastic differential equation 

(6.1) dZ.^t = o,{Zi^t) d^M + ^Zi^t dt + F^{Z^) dt 1 < i < iV, 

where the <T(zj) are the diffusion m x m-matrices, the ^j^t are independent 
standard Wiener processes valued in M™, ^ is an m x m-matrix and the 
jpN . -^m _^ jgm force fields acting on each particle. Because of 

indistinguishability we assume 

ry\ TpN I ^ ..N-1 



with Zj^ := {zi,..., Zi+i, ...,zn) and : M'" x T'(M'") ^ M"*. (Note 
that here and below the latin letters "i, j, ..." label the particles, whereas 
the greek letters "a, /3, ..." label the coordinates). 

We assume that is uniformly bounded and Lipschitz in both variables 
(when endowing V(M"^) with a distance inherited from a negative Sobolev 
norm). More precisely, we assume that for any k > m/2 there exists CF,k > 
such that for any z,2e M™, /, / G P(M™) 

(6.2) ViVGN*, \F''iz,f)-F^iz,f)\<CF,k[\z-z\ + \\f-f\\H 

It is also natural for the limit to exist to assume that there exists a function 
F : M"* X P(M™) such that F^ F, in the sense that there is a 

constant Cp^nm > such that 

(6.3) VAf G N*, Vz G M™, V/gP(M™), \F^ {z, f) - F{z, f)\ < 



N 



F- [Z, ^--^ j = ^ F- [z,, Z,) , F- [z,, Z,) ■■=j^Y,^{z,- z,) 



A simple example which satisfies these assumptions is 
(6.4) 

( ..N-l\ _ ^ TPN ^\ CnAf / V \ ._ 1 , 
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for a smooth vector field ^ : ^ M*", so that F{z, /) = * f){z). 

Under the smoothness assumptions (16.2^ on the A^-particle force fields, for 
any A'' > 1 there exists a Markov process {Zj^)t>o which solves the system 
of stochastic differential equations (|6.ip , see |371 130] . 

The time-dependent law //^ of the process Z^^ satisfies the following linear 
master equation corresponding to (j6.ip . given in dual form by 

(6.5) V(^GP(M'-), 9i(//^,vp) = (//^,G^<^> 

where is defined by 

N N 

i=l i=l 

N 

i=l 

The nonnegative diffusion matrix A, the gradient V, and the Hessian matrix 
V? associated to the variable Zj = (-Zj^i, . . . , ^i^m) ^ corresponding to the 
i-th particle are given by 

^ 7=1 



and 



We also introduce the nonlinear mean field McKean-Vlasov equation on 
): 



(6.6) ^/ = g(/,), /|,=o = /i„ in P(M-), 
with 

m mm 

Qif)= E dlf,{A^,pf)-Y,da[i^zU]-Yda{F^{zJ)f). 

a,/3=l 0=1 a=l 

There is an important literature on this class of nonlinear partial differ- 
ential equations. See fore example |28) and [8j where more details and more 
references can be found. 

In the sequel, we make the following strong structure, smoothness and 
boundedness assumptions on the coefficients: 

(6.7) {A = 0, k:=0) or {A > kU, k > 0, A G Ty^'~(M™)), 
as well as 



(6.8) Vz G M'", V/ G P(R"), F{zJ)= '^{z-z)f{dz) 
where ^ G F|''(M™) for some k eN, k > m/2 + 3. 
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When the diffusion matrix A = is zero, m = 2d, z = {x,v) £ M^^, our 
assumptions cover the case of the mean-field Vlasov equation. Indeed the 
classical Vlasov equation reads 

dtf + v-V.J + {V.,4^*p[f])-VJ = 0, f = f{t,x,v), x,veR'', 



with 



P[f]{t,x) = / f{t,x,v)dv, 



and it falls into our structural assumptions with z = {x, v) G M'^ x 
^(z) = -^{x) = {0,V^tlj{x)) with V^V' is i^^^+^+o ^^^^ 

Then F = {Fx,Fjj) defined by (|6.8p is given by 

F^{x, v) = 0, F^{x, v) = V^^P * p[f] 

for the limiting system, and, with X G (M'^)^ and V G (W^)'^ , we have for 
the A^-particle system ^ defined as above and F^ = {F^ ,Fy) given by 



F# = 0, (F^), = 1^V,V'(X,-X,), i = l,...,N. 

Observe in particular that it does not allow for the Coulomb or Newton 
interactions in this Vlasov setting due to the smoothness assumption on ip. 

6.2. Statement of the result. Our main result in the section is a quanti- 
tative propagation of chaos result for the class of equations described above. 
We state two separate results respectively for the McKean- Vlasov case (pos- 
sibly non-zero diffusion matrix) and the Vlasov case (zero diffusion matrix). 

Theorem 6.1 (The McKean- Vlasov equation). Consider an initial dis- 
tribution /in G Vq{W^), q > 2, the hierarchy of N -particle distributions 
fl^ = S{^{f^^) following (j6.5p and the nonlinear evolution ft = ^^^(/in) 
following (j6.6p . Assume that (|6.7p and (j6.8p hold. 

Then there is k £ N and a constant C > and, for any T > 0, there are 
constants Ct, Ct > such that for any 

^ = ipi®---® ipt^T^^ J-:=/7|(M'")nLip(R"^), \\ipj\\T<l, 

we have for N > 2i: 
(6.9) 

sup|((5f (/^) - (5f^(/i„))^"^) ,^)\ <C^^ + Cr^^ + CTinT{^^)- 

As a consequence of (j6.9p and Lemma \4-()\ this implies the propagation 
of chaos with rate e{N) < C{£, T, /;„) for any initial data /;„ G 



36 



S. MISCHLER, C. MOUHOT, AND B. WENNBERG 



Now we consider the case of the Vlasov equation. As will be clear from the 
proof, when ^ = and (0) = 0, the error e(A^) = vanishes in assumption 
(A3). This leads to the following improved result. 

Theorem 6.2 (The Vlasov equation). Suppose, in addition to the assump- 
tions for Theorem \6.1\ that ^ = and '^(0) = 0. Then there is a constant 
C > and, for any T > 0, a constant Ct > such that for any ip G Lip(M^'") 
and any N > I: 



(6.10) 

sup 

[o,r] 

(observe the replacement ofW2 by Wi in the last term) which in turn implies 
[0,T] iV \ / 



^C^Ct 07 (/in) 



N N 



Remark 6.3. Note that the coupling method introduced in [37\ leads to a rate 
of chaoticity of order 0{1/^/N) for the normalized Wasserstein distance W2 
between the law of and the tensor product f®^ . This is better than 
our estimate, which is limited by the estimate in Lemma \4-6\ However, the 
coupling method is usually limited to the quadratic interaction given by (16. 4p . 



6.3. Proof of Theorem 16.11 As in the proof of Theorem lS.ll we prove that 
Theorem 16.11 is a consequence of Theorem 12. H by verifying that assumptions 
(A1)-(A2)-(A3)-(A4)-(A5) hold. However, in the present model we can- 
not, as in Section \5\ use the total variation norm for the key consistency 
estimate (A3) and differential stability estimate (A4). The reason is that 
qN^N^ involves derivatives of Z 1— t- ^{(j)^), and hence of Z 1— )■ /i^ which 
is not differentiable from M""^ to P(M™) when P(M'") is endowed with the 
total variation norm. We therefore make the following choice of functional 
spaces: E := with 

Ti = H^'{R"'), si>f + 2 

and the weight rng-^^{z) = rag^{z) = 1, and = Lip(M™) and Vg^{E) := 
7^2(1^"^) endowed with the quadratic MKW distance W2- 
Proof of assumption (Al). The symmetry assumption is a consequence 
of the fact that first (|6.5p is well posed for any f^ G Vsyra{W^^) so that 
//^ is a probability measure for any t > 0, and second the generator 
commutes with the permutations. 

Proof of assumption (A3). We claim that for any si > m/2 + 1 there 
exists a constant Cs^ such that for all ^ G C^'^(T'g^ (£^), M) 

(6.11) ||G^(<I>o^f)-(Q(^f),D$[/.f])||^^(^,)<^||cI>|b.,.(p,^(,,),M), 
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Proof of (j6.1ip . First, the map 
is with 

Take $ G Cl^iVg^iE)). Then the map 

is Cfe^. Indeed, denoting (/> = (/.^(•) = D^[^i^] G (FI2'HIK"))' = ^2'(IK"), 
we can write: 



and both 5^- ^f^^. and ^ ^"^z^ belong to ff_2^(M"*) thanks to the condition 
si > m/2 + 2. 

As a consequence, we compute 

(gV$) {z) = g'' ^^l^z) 

N 
1=1 

N N 



/l(Z)+/2(Z) 



with 



i=l o,/3=l 

N m N m 



A 

=1 a,^=l i=l a=l 
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and 



N m 



1=1 a,/3=l 
N m 

i=l a=l 



+ 



On the one hand, using that 



2 C 

< T7 sup IIOzJIj^-sWIRm) < 



as weh as (|6.2p and (|6.3p . we deduce that 

2 

\I2{Z)\ < Nj^\\A\\^\\DH\\^\\d,5\\l^s, 



H 



On the other hand, we recognize 

/ m 

h{z) = (f^^, ^'^A^z 



C(j) 

7^' 



\ Q=l / \ a=l / 

(Q(^f ),<^z> = ),D1>(/if )> = (7r^G-<I>) (Z), 



thanks to the calculation of the limit dual generator made in Subsection 14.11 

□ 



Proof of assumption (A4). We need here to perform a seconcZ-order 
expansion of the limit semigroup. 
We consider 

• for any two given initial data /;„, g-^^ G V{W^') the corresponding solu- 
tions fi and gt to the nonlinear McKean-Vlasov (or Vlasov) equation 

dSS]), 

• for any given initial data € V{W^^ the solution ht to the following 
equation, which is the linearization around ft'- 

(6.12) 

dth = V2 : (A/i) - V • aSTz) h)-V -[hi^ *f) + f{^ *h)], ht=o = /Jin, 



(6.13) 
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rt the solution to the fohowing second variation equation around ft 



dtr = : (^r) - V • {{STz) r) - V • [r * /) + / * r) 



1 

-2^ 



hi^ *h) 



1 

2^ 



hi'i/ *h) 



r\t=Q = Tin = 



for two solutions /i, h of the first variation equation. 
Then we shall prove the following a priori estimates. 

Lemma 6.4. For any si G N, si > m/2 + 1, £ G {1, 2, 3} and for any T > 0, 

there exists Ct such that 



.14) sup \\gt - /tll^-^imm-) < Ct hin - finW^-^l 



[0,T] 



3.15) sup ll^tll^-n(Km) < Ct ||/ii„||^-n(jj™), 

3.16) sup < Ct ||/iin||^-^i 



[0,T] 



II/). II 

m) ll"in||^-n(]j,„-), 



anc? if/ien = hm = g-m — fm we have 

(6.17) sup \\gt -ft- ^ Ct \\gin - fh 



[0,T] 

sup 

[0,T] 



|2 

'/i'l2l(IR™)' 



3.18) sup \\gt - ft-ht- n||^-{.i+2)^j^„^ < Ct \\gin - fh 



|3 

'_H"l2i(R™)' 



This shows that the nonlinear semigroup S^^ associated to the nonlinear 



McKean-Vlasov equation 116. 6\) is C'^'^{Vg^{E),Vg2{E))■ 



Proof of Lemma \6.4\ The proof is carried out in several steps. 
Step 1. We will several times consider the equation 

(6.19) dtQt = V2 : (AG) - V • ((=^z)G) - V • (m G + U2 (^ * Qt)) 

with given initial data and with an M™- valued function ui and an M- valued 
measure U2 to be specified (chosen in order to "match" equations (16. 6p . (I6.12P 
and (|6T3]) ). We claim that for any A: G N, A: > m/2 + 1, £ G {1,2,3} and 
any T > 0, 



(6.20) Mt G [0,r], IICtllHI^fe(Km) < IICin||j:^-fc(Km) 6 

with 



Ck{'^ ,ui,U2) ■= C{k) sup 

te[o,T] 



40 



S. MISCHLER, C. MOUHOT, AND B. WENNBERG 



We argue by duality and we consider a smooth solution 9 to the following 
linear equation (which is the dual equation of ()6.19p ) 

' dte = L\e + Lie, 
(6.21) \ Lie := A-.v'^e + {,^z) -ve, 

^ Lie ■.= ui-ve + ^ * {u2 V9) , 

with "^{x) := '^{-x). 

For a given multi-index G N™ with {I'l = k' < k, we compute 

A / \d^e\'Uzrdz 

at J Km 

= / {d''Lie)d''e{zfUz+ [ {d''Lie)d''e{zfUz=:Ci + C2. 

By integrations by parts, we get 

A<-K^ |Vrenz)2^dz + C (||=^||oo+ P|lH/fe'+2,oo(R™)) ll^llJ^*'(Rm) 

and (using Sobolev embedding inequalities for the last term) 

-^2 < C (^ll'"l|lvi/'='.°°(R'") + ll'"2||Ty(IR™) ll^ll_fffe'(R™.)) ll^*ll/i'*{R'") 
which shows that 



Vt G [0, T], 116**11 wfefR"!) ^ 



) ^ \\"ia\\H^{Rm) e 

Denoting by Ut the linear semigroup associated to (j6.19p . the associated 
dual semigroup U* is generated by (j6.2ip . As a consequence, for any E 
H^{W^), we have 

(Ci)^in) = {CimUf Oin) < || Cin || j;^-*: (jjm) I l^t ^in 1 1 //|(Rm ■) 

< Cfee'^,«l,n2)T 11^. II |m II 

— W^mWjJ-k^^n^-^ IrinllHfe(Rm), 

which concludes the proof of the claim (|6.20p . 

Step 2. Proof of (16. J4|) . The equation satisfied by the difference dt = gt — ft 



Cfc{'2<,«i,«2)T 



IS 



(6.22) 

r dtdt = V2 : (^ di) - V • {i^z)dt) -V-{dti^*f)+gi^* dt)) . 

[ d|t=0 = din = ^in — fin, 

which fits in the form (j6.19p with ui := ^ * / and U2 = g. Now, since 



V\^*f) 
we conclude that 



L°°(R'") 



iV'"^)*f 



< 



L°°(R'") 



L°°{R'") 
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and that (I6.14P holds. Proceeding in the same way for the function h we end 
up with 



(6.23) 



sup \\ht\\jj-si^^m) < Ct WhinW^-jif^r^Y 



for any si G N, si > m/2 + 1. 

Step 3. Inequalities for products and convolutions in Sobolev spaces. We 
define the weighted L°°-based Sobolev spaces as usual: 



0<k'<k 



We have the three following inequalities on functions S, ip in the appro- 
priate spaces: 

\\S *'^\\l°^^{W^) < ||^||//fe(]Rm) ll'S'llj^-fc(Km) 

and more generally 



15*^1 



< \\^\ 



and finally 



\H- 



< Chj \\s 



k,i \\'~>\\h- 



for any k,i £ N, and some constant C^^i > 0. The proofs are elementary 
and we omit them for the sake of conciseness. 

Step 4. Proof of id. 17]) . Let cot '■= gt — ft — ht = dt — ht, which satisfies the 
equation 



in 



0, 



(6.24) 
with 

r Luj :=V^ : {Auj) - V ■ i{^z)ujt) - V ■ {u {"^ * f) + f {"^ * u)) , 

1 St = V-(da^*dt)). 

Denoting by Qs,t'w the unique solution of the linear, non-autonomous 
equation 

dtwt = Lwt, Ws = w, 
the Duhamel formula for equation (|6.24p yields 

rt 



J 

Jo 
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Therefore we obtain, using (I6.20p and the estimates estabhshed in the Step 3, 
that for any t E [0, T] 



\^t\ 



H- 







< Ct 



|V(d, (^*d. 



ds 



iVd 



-2 V"^ ) 
_ t 



* ds||„,fe,oo,^™^ ds 







1^ * (Vd, 



< CtM\ 



-(k-i)f 



ds, 



which together with (|6.14p for the control of the norms of dj implies (|6.17p . 
Step 5. Proof of (j6.16p and i6.18\) . The second variation r satisfies the 
equation 

dtr = Lr + Rt, r|t=o = 0, 

with L as above and 



R. 



1. 



t ■- 



1. 



V ■ (htC^ *ht)) +-V ■ (ht i'^*ht) 



'H 



-(k-i). 



We proceed as in Step 4, taking advantage of the bound (16.23p . and we obtain 
(6.25) sup Iktllj^-fe(ffim) < Ct ||^in||^-{fc-i)(]g„) WK, 

which is nothing but (j6.16p . 

Finally we introduce ipt '■= gt — ft — ht — rt = dt — ht — rt = ut — rt, with 
the initial data /lin = /lin = din — fin- It satisfies the equation 



(6.26) 

Therefore, we deduce 

yte[o,T], I 



H- 



< 



t ds 



< Ct 



f 

Jo 



-(k-i). 



+ lid. 



-(fc-i). 



ds, 



which together with (f6mp - ([6T5]) and (f6TT71) implies (f6TT8]l . 



□ 



Proof of (A2). The first property (A2)-(i) is a consequence of (j6.14p in 
Lemma \6A\ we have 



(6.27) 



IIQ(/l)llHl2*^(Rm) < C^^^l 
||Q(/2)||j^-2(Km) < Co^^l 

\\Q{f2) - Q(/i)IIj,-^^(m™) < C.^,2 11/2 - /l||^/'.( 
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We write 

Q{f) = Qi(.f) + Q2{f) 

with 

Qi(/) = V^:{Af)-V- {{^z)f), Q2U) = -V • ((^ * /) /). 

The Unear term Qi satisfies the first part of (|6.27p (first equation) by direct 
inspection combined with the use of Sobolev embeddings. In order to see 
that Qi also satisfies the second part of (j6.27p (on the difference), we write 



||Ql(/2) - Ql(/l)||j:^-fe(IRm) < \\Mw'^-'^'^{9."^)\\h- h\\jj-(M-'i)(^^rn) 

+ ||=^||oO 11/2 - /l|l^-{fc-l)(]i5m)' 



and we conclude by using interpolation and Sobolev embeddings (noticing 
that k — 2 — 1/2 > m/2 allows for the Sobolev embedding in the first term). 

Concerning the quadratic term Q2, using estimates proved in Step 3 of 
the proof of (A4), on the one hand we get 

ll<52(/)||j:^-fe(IRm) < ll<32(/)||H-fe(R™) < 11'^ * /lll^fe-i.'^CR™) ll/llH-(fe-i)(IR™) 

< Ck ||'^||j^2(fe-l)(IRm) ||/||^_(fe_l)(.]g™j < Ck ||'^||j^2(fe-l)(]gm) 

where we have used 7^(1^"^) C H~^''~^^^'^{W^) with continuous embedding. 
On the other hand, we have with d := /2 — /i 

||Q2(/2) - <52(/l)||j:^-fc(Km) < Wi'^ * d) (•)"^||tyfc-l,oo(]Rm) 1 1 /2 1 1 J^- (fe- 1) (Rm ) 

+ W'^ * /l|lw'fe-i.°°{K'") ll^ll/f-('=-i)(Mm)- 
In order to estimate the first term in the above inequality, we remark that 

|(a"^ * d)(z)| < wd'^^iz - .) ||d ir'WH-HR-) 

uniformly for any z S M™. All together, we have for the quadratic term 

||<32(/2) - <52(/l)||j:^-fc(Igm) < 11/2 - /lllHI^'=-l)(IRm)' 

and we conclude the proof of (A2)-(ii) by using interpolation and Sobolev 
embeddings again. 

Proof of (A5). We use the well known following estimate (see |37)): for 
any q > 1, fin, Sin S Vq{M.'^) and T > there exists Ct such that 



supl^,(5f^(/in),5f^(gi„)) < CTW,iU,9in), 
t>0 

that we use with q = 2. Alternatively, estimate (j6.14p precisely says that 
assumption (A5) holds in Vg^{E). 
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7. Inelastic collisions with thermal bath 

7.1. The model. In this section we assume that E = 'R'^, d > 1, and we are 
interested in the following Boltzmann equation for diffusively excited granular 
media on the distribution f{t,v) > 0, v £ M.'^ oi particles: 

(7.1) ^ = Q{f,), /|,=o = /in in V(R^), 

with 

Q(/) = Qa(/,/) + ^A/ 

for some > 0, and where the quadratic Boltzmann collision kernel Qa is 

defined by the following dual formulation 

(7.2) 

{QM, /), ■■= ! Hcos 9) {<P{wl) - ct>{w2)) da fidwi) f{dw2) 



^ WI+W2 U* 



1 — a\ / 1 + a\ 

[Wi - W2) + \W2 - wi\ a, 



for any (p G Co(M°), / e 'P(M ), and with cos9 = a ■ {w2 — Wi)/\w2 — wi\ 
and similarly as in equations (jS.ip and (j5.6p 



but with 



for some a G (0, 1). (Note that the case a = 1 and u = {) would correspond 
to the elastic Boltzmann kernel considered in Section O) This corresponds 
to a situation where particles lose energy when they collide. We refer to 
[21 [7] for a physical motivation to these equations. The mathematical theory 
is treated in e.g. j3l[2l|4|, where, for example, it is proven that this equation 
generates a nonlinear semigroup S^^fin '■= ft for any /in S Vq{M'^), q > 2. 
Notice that unlike the classical Boltzmann equation the kinetic energy is not 
conserved. For the sake of simplicity we make the normalization assumptions 



Ili(S''-i) = / b{ai)da = l, u = I. 



One of these quantities, say the first, can be set to one just by a rescaling 
of time but the two cannot be changed independently. However, the result 
would be the same for any value of u. Note that due to the normalization 
||6||2,i(gd-i) = 1 and the fact that / G V{M.'^), the bilinear operator Qa splits 
into a quadratic part and a linear part 

Q{f) = Qt{fJ)-f + ^f, 

where is defined through the positive part of the expression (j7.2p . 

We now want to introduce a A^-particle system associated to the above 
Boltzmann equation for diffusively excited granular media by mimicking the 
Kac's construction. We consider the velocities process (V/^) with values in 
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W^^ , of mixed jump and diffusion nature, defined through the stochastic 
differential equations 



N 

Here is a M"^^ valued standard Brownian motions, AA^(ds, dcr, i, j, dz) 
is a Poisson measure on [0, oo) x S*^"^ x {1, . . . , N^'^ x R_|_ with intensity 



/ / E ^^^U^l ) l.<6(..n. , (V- )) -^^(^^^ da, j, dz) + . 



1 ^ 



independent of and the two functions ^i^j^a '■ M'^^ — t- R'^^ and Uij : 

^dN _^ §"^"1 are a.e. defined through the following expressions: for any 



V = {vi,...,vn) GR'^^ we set 



Ui 



and 
where 

and, as in equation (15. ip . 

(7.3) < = f + f. .- = ^-1. 

but here with 

1 — q\ / 1 + a 



Wij = Vi + Vj, = [ ] 



N 
t 



The associated forward Kolmogorov equation on the probability law /, 
of (Vf ) in R''^ reads 

(7.4) 5i(/,^,(^) = (/,^,G^(^) 

with generator = Gi +G2 , where Gi is associated to an inelastic Boltz- 
mann collision process whose collision kernel only depends on the deviation 
angle as in Section [5] 

1 ^ r 

(7.5) {G^^m = ^ E / ^(^o^^^J-) ['^(^iP - '^(^)] 

with cos 9ij = a ■ {vj — Vi)/\vj — Vi\ and V-* defined in (|7.3p . and is the 
generator associated to the Brownian motion 

N 

(7.6) {G^^){V) = Y,^^^, 

i=l 
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where Vi := . . . , Vi^d) and Aj denotes the Laplacian in M"' associated to 
the i-th particle: 

d 

a=l 

It is classical to prove that V/^ is a Feller process and we refer to the textbooks 
[13| l33:j where the theory is set up with full details (one can also refer to 
[37 | 130 1 [T^ where similar processes are considered). 



7.2. Statement of the result. The main result in this section is a quanti- 
tative estimate of propagation of chaos for the mixed collision and diffusion 
model introduced above. 



Theorem 7.1. Consider an initial distribution fin E Pq(M ), q > 2, the 
hierarchy of N -particle distributions //^ = 5'j^(/^^) following the evolution 
()7.4p . and the nonlinear semigroup ft = 5*^^ (/in) following the evolution 

m- 

Then there is a constant C > and, for any T > 0, there are constants 
Ct, Ct £ (0, oo) only depending T G (0, oo) such that for any 

we have for N > 2i: 
(7.7) 



sup 

[o,r] 



As a consequence of (j7.7p and Lemma \4-6\ this shows the quantitative prop- 
agation of chaos with rate e{N) < C{i,T, f^^-^) for any initial data 
/i„GPrf+5(M^)- 

We are not aware of any result of propatation of chaos in this setting. 
A conceivable alternative approach would be to use the general nonlinear 
martingale approach, but that would most likely not provide any quantitative 
rate of propagation of chaos. The techniques developed recently in |15| for 
the elastic Kac equation without cut-off is yet another alternative technique 
that could be tried on this model, but we have not made any attempts in this 
direction, and, would it work, it is not clear as to what kind of convergence 
rate one could hope to achieve. 



7.3. Proof of Theorem 17. IL We shall prove that Theorem 17.11 is a conse- 
quence of Theorem 1 2 . 1 1 bv proving that the assumptions (A1)-(A2)-(A3)(A4)- 
(A5). We consider the phase space E = M."^ and the following choice of 
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si := 3, 
S2 := 3si 



9, 



where the Fourier based space % ^(M ) and the norms | • 1^ are defined in 
example 13.71 and the corresponding spaces Vg^ [E) and Vg^ i^) (without 
weight). We finally define ^"3 = Lip(R'^) and Vg^iE) := P2(M'^) endowed 
with the quadratic MKW distance W2- 

Proof of (Al). The well-posedness of equation ()7.4p - ()7.5p is a variation 
on the well-posedness result for equation (j7.ip as obtained in |[2| |4j . We also 
refer to \13\ [571 EHl US] for a proof of the fact that t is a Feller process. 

Proof of (A2). First we prove (A2)-(i), and more precisely we prove that 

sreC^'HVgAE),Vg,m, 
which is a consequence of the following result: 



Lemma 7.2. For any fin, gin G V{W^) and any final time T > 0, the asso- 
ciated solutions ft and gt to the diffusive inelastic Boltzmann equation (|7.ip 
satisfy for any s > 

2T 



(7.8) 



sup l/i 

te[o,T] 



9t\,<e'' l/i. 



5ir 



Proof of Lemma We recall Bobylev's identity for Maxwellian inelastic 
collision kernel (see for instance [2]) 

^{QW\g)){0 = Qt{F,Gm=:\ I bU-i) [F+ G- + F- G+]da, 
with F = f, G = g, F^ = F(e±), G± = G(C±) and 



(e-iek). 



4 - ' 4 ' - 4 

Denoting D = g — f , S = g + f , the following equation holds 
dtD= I h{a-i 



Using that ||5||oo < 2 and then < |^|, we deduce in distributional sense 



(7.9) 



D+ S- D- 5^ 
h — 



da-D-\C\^D. 



d \D\ 



< 




da 



< 2 sup 



from which we conclude that (17.811 holds. 



□ 
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Next we prove (A2)-(ii), as a consequence of the following result: 
Lemma 7.3. For any f,g^ V{W^) and s > 0, we have 

(7.10) \Qa{fJ)\s<2 

and 

(7.11) \QM + 9J-9)\s<^\f-9\s. 
Moreover for any s > 2 there exists (5 E (0, 1) such that 

(7.12) IA/-A5I, <2|/-5|f. 

Proof of Lemma \7.SX We prove the second inequalities (IT.lip . We write in 
Fourier: 

J'{QM' + 9,f-9))=Qc.{D,S) 

= \ ! K^-i) {S{e)Dir) + S{r)D{t)-2DiC)) da 
^ J§<i-i 

where Qa is the Fourier transform of the symmetric version of the collision 
operator Qa, which yields 



Qa{D,S) 



<Ti + r2 + T3, 



with 



Ti 



2(0^ 



< 



bia-0 



\sie)\ mm ic 



da < \D\ 



2 iO' 

Similar estimates hold for the two other terms T2 and Ts- The proof 
of the first inequality (IT.lOp is similar (and simpler): we use the Fourier 
representation of Qaif,f) and the bound < 1. We finally prove the 

last inequality. We compute 

<5\ 



|A/ -Agl= sup leP^^^ < sup ( |F - G\'-' 



\F-G\ 
with (5 := (s — 2)/s. 

Proof of (A3). We claim that for any Si > 3 there exists Ci G 



□ 
such 



that for aU $ G C'^'\Vg^iE),R) 
(7.13) 

||G^(cI>o;.f)-(Q(^f,;.f),I)cI>[;.f])ll^^(^,^ _ ^ 

which is (A3) with A; = 2, r/ = 1 and e{N) = Ci N'^ . 

We begin with a technical lemma which shows that the norm | • |s is 
well-adapted for obtaining differentiability of the empirical measures. It is 
worth emphasizing that the choice of si = 3 (in fact we only need si > 2 
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by modifying slightly the arguments) comes from the requirement that the 
function V ^ he C^. 



Lemma 7.4. The map 



t,dN 



Vg.iE), fi^ isC^'^ and 



Proof of Lemma \7.4\ For u , E M*^, we have 



|5„ - 5^ 



sup 



(0^ 



<\v — w\ sup ■ 



< \v — w\ sup — — < \v — w\ 



which shows that w i— )• (5^ is C^'^. For the sake of simplicity we present the 
proof of differentiability when d = 1, the case d> 1 being similar. For f G M 
and /i G M*, we have 



< sup 



from which we deduce that v ^ 5^ is C^'^ . Similarly we can go to second 
order: 



sup ■ 



<SUp^<|/.|3, 



and we easily conclude that v ^ 5^ is C^'^. When the dimension d is greater 
than 1, one can perform the same argument for the partial derivatives of the 
Dirac mass. □ 

We come back to the proof of (I7J3I1 . Take ^> G C'^^^{Vg^{E),W) and 
compute separately the contributions of Gf ^ i = 1,2. Proceeding as in the 
proof of (A3) in Theorem 15.11 we have 



o ) = (Q, (^f , ) , Z)<I>(^f )> + I,{V) 



with 



1 ^ r 

\h{V)\ < ^ E / Kcos(%))||f I 



I] 



dfj 
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since for any i ^ j 



My.* 



Sl 



< 



1 

1 



+ 



3 ISl 



4 



On the other hand, as in the proof of assumption (A3) in Sectional the 
map W'-^ — > M, y i-> is C^'^ thanks to Lemma 17.41 and denoting 

(pz = -D$[//f] e we compute 

AT 



i=l 



+ 



l^llc2.i(Pg^(£;),i 
iV 



□ 



We conclude the proof by combining the previous estimates. 

Proof of (A4). For /in,5in £ 7^(M'^), we define the associated solutions ft 
and gt to the nonlinear Boltzmann equation; we define ht := [/in] (^in — 
/in) the solution of the linearized Boltzmann equation around ft] and we 
define rt the solution to the "second variation" equation around ft- More 
precisely, we define 

r dtft = QMuft) + ^fu f\t=0 = fin 



dtot = Qaigt,gt) + ^gt, g\t=o = gin 

dth = Qaift, hi) + Qa{ht, ft) + A ht, h, 
^ dtn = Qaift, n) + Qain, ft) + A rt + \Qa{ht, h) + \Qa{ht, ht), r\t=Q 



i|t=0 — ^iii) 



where in the last equation (second-order variation) ht and ht are two solutions 
to the third equation (first-order variation). 

We then define when h^^ = h\a = 9in — /in the following error terms 

r dt := gt - ft 

^t := gt-ft-ht = S^Hg-n.) - 5f^(/i„) - /:r[/i„](5in - /in) 

. ipt ■= gt- ft - ht- n. 

Lemma 7.5. Fix s > and T £ (0, oo). There exists Ct such that for any 
fin, gin £ V{M.'^), the following estimates hold 

(7.14) VtG[0,T], <Cr |/iin|,, 



(7.15) 



VtG[0,T], \rt\2,<CT\h 



in 1 5 |^in|s; 
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and when hin = ^in = ^in — fin we have furthermore 

(7.16) VtG[0,r], |a;t|2, <Ct |/in-5i„lL 

(7.17) VtG[0,r], |Vt|3s<C^T |/in-5in|'. 

This proves that Sf^ E C'^^^{Vg^{E),Pg^[E)). 

Proof of Lemma I7.5i We skip the proof of (|7.14p since it is similar to the 
proof of (|7.8p . We then deal with each term successively. We work in Fourier 
variable and we introduce the notations F = f, D = d, H = h, H = (h), 
Q = (2}, R = f and ^ = ip. 

Step 1. The evolution equation satisfied by is 

(7.18) dtn = Q^{n, F) + Q^{F, n) - i^p n - q^{d, d). 

We deduce in distributional sense 

d m)\ 



dt (0^ 



< Ti+Ti, 



where 
Ti :-- 



sup 



2s 



n{c+)F{c 



+ 



< sup / b [a ■ ^ 

sd-l 



\2s 



2s 



+ 



m~)\ (t 



+ \m\ 



da 



< c sup lif + |n(o)| sup I'^K*! 



for some constant C > 0, and 



1 

- sup 



2 ^giRdy§d-i (0 



2^ iD(e+)i?(r) + ^(r)^(r)| da 



< 



sup 



6 a-C 



iz)(e+)i |i)(r 



using the estimates (j7.8p . We then conclude thanks to a Gronwall lemma. 

Step 2. The evolution equation satisfied by R is 
(7.19) 



1 A 



1 



a^i? = Q,(F, i?)+Q„(i?, F)- lel' R+-QaiH, H) + -Q^iH, H), i?u=o = 
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Equation (I7.19P being similar to equation (j7.18p , with the same computations 
as in Step 1 we deduce that (j7.15p holds. 

Step 3. Choosing now = = din, the equation satisfied by ^ is 

Observe that by conservation of mass ^'j(O) = for all times, and with these 
choices of initial data Qai^,H) = Q+{i},H) and Qa{D,n) = Q+{D,n). 
Then we perform similar computations as in Step 1, and we deduce in dis- 
tributional sense 

d 1^(01 



dt (0 



3s 



where 



Ti := sup ^-3^ < C sup ---3^^ 



_ \Q+{^,H)\ ( ( \H{i)\ 



\Qt(D.a)\ \D{m\ I \m\\ 

75 ■= -P ^ip^ < 2 sup ^ j ™p ^ j 



Finally we then conclude the proof of (j7.17p using the already established 
estimates (|7.8p . (j7.14p . (|7.16p . and the Gronwall lemma. □ 

Proof of (A5). We use the following result proved in [3] (see also [2] for 
a similar result) 

SUpVF2(5f''/i„,Sf^5in) < Ty2(/in,5in), 

which concludes the proof. 
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